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ABSTRACT 

The characterization of the wind speed vertical distribution V(h) is fundamental for 
an astronomical site for many different reasons: (1) the wind speed shear contributes 
to trigger optical turbulence in the whole troposphere, (2) a few of the astroclimatic 
parameters such as the wavefront coherence time (to) depends directly on V(h), (3) the 
equivalent velocity Vb, controlling the frequency at which the adaptive optics systems 
have to run to work properly, depends on the vertical distribution of the wind speed 
and optical turbulence. Also, a too strong wind speed near the ground can introduce 
vibrations in the telescope structures. The wind speed at a precise pressure (200 hPa) 
has frequently been used to retrieve indications concerning the tq and the frequency 
limits imposed to all instrumentation based on adaptive optics systems, but more 
recently it has been proved that V200 (wind speed at 200 hPa) alone is not sufficient to 
provide exhaustive elements concerning this topic and that the vertical distribution of 
the wind speed is necessary. In this paper a complete characterization of the vertical 
distribution of wind speed strength is done above Mt. Graham (Arizona, US), site of 
the Large Binocular Telescope. We provide a climatological study extended over 10 
years using the operational analyses from the European Centre for Medium-Range 
Weather Forecasts (ECMWF), we prove that this is representative of the wind speed 
vertical distribution at Mt. Graham with exception of the boundary layer and we 
prove that a mesoscale model can provide reliable nightly estimates of V(h) above 
this astronomical site from the ground up to the top of the atmosphere (~ 20 km). 
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1 INTRODUCTION 

The characterization of the wind speed above an astronom- 
ical site is extremely important for several different reasons. 
Firstly, because the wind speed is strictly correlated to the 
optical turbulence: strong wind speed and sharp wind speed 
gradients are indicators of a turbulent atmosphere, which in 
combination with a stable stratification of the atmosphere (a 
positive gradient of the potential temperature) creates op- 
tical turbulence that limits the resolution of the telescopes. 
Secondly, because the stronger the wind speed, the higher 
is the speed at which the turbulence layers cross the pupil 
of the telescope and the higher is the frequency at which 
adaptive optics systems are forced to work to correct the 
turbulence perturbations above the wavefront. Thirdly, be- 
cause a too strong wind speed near the ground can introduce 
vibrations of the telescope structures. 

At mid-latitudes, the wind speed varies with height 
showing a maximum at the jet stream level, usually 10-12 
km (or ~200 hPa) above sea level. The strength of the wind 
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speed at high altitudes varies according to season, with the 
strongest values during (the local) winter and early spring 
IIMasciadri fc G arfias 20 0ll;|Carrasco. Avila fc Carraminana 



2005 : iGarci'a-Lorenzo et all [20051: iMasciadri fc Egnerl | 200 
Eerier, MasciadrTfc McKenn<Tl2007: Bounhir ct al. 200 



Masciadri et al.ll2010l ). It has been observed that seasonal 



variations at the jet stream level are more pronounced at 
higher l atitudes than at low latitudes in proximity of the 
equator (jCarrasco et a!ll2005T ) because of the general circu- 
lation of the wind speed at synoptic scale. 

At extreme latitudes, i.e. in proximity of the poles, 
the wind speed is characterize d by a completely different 
feature in the free atm osphere IjGeissler fc Masciadrill2006l : 
ISadibekova et all 120061 ). It is rather weak up to 10 km and 
then increases monotonically above this height with a rate 
that increases with t he distance of the s ite from the centre 
of the polar vortex (Hagelin ct al. 2008). This assumption 
has been confirmed bv lLascaux et all ([2009). 

Apart from the fact that a windy atmosphere more eas- 
ily triggers optical turbulence, the astroclimatic parameter 
that directly depends not only on the optical turbulence but 
also on the wind speed is the wavefront coherence time: 



2 S. Hagelin et al. 



0.057A 6/5 (' J V{h) 5, ' A C 2 N {h)dh\ 



-3/5 





also equivalent to: 
ro 



r = 0.31 



Vo 



(1) 



(2) 



where ro is the Fried's parameter and Vo is the equivalent 
velocity: 

3/5 



/ 'Jv(h) 5/3 c 2 N (h)dh S) 



V = 



\ 



JC%dh 



(3) 



/ 



to depends, therefore, on the vertical distribution of 
the wind speed V(h) and the optical turbulence C% (h). 
To simplify the calculation of the wavefront coherence time 
Sarazin & Tokovinin (2002) proposed a method to calcu- 
late to using only the wind speed estimated at 200 hPa in- 
stead of the whole profile using an empirical relationship 
between Vo and V200: firstly they used Vo — 0.4 * V200 
and, in a second time, the expression has been modified in 
Vp = m.ax(V orourl d,0A * V200). However, a few years ago 
IjMasciadri fc Egnerl 2006), in a study done above San Pe- 
dro Martir, it has been observed that the value of the con- 
stant 0.4 was not universal. The authors found the value 
0.56 above San Pedro Martir. However, they also proved 
that the relative error intro duced in to with the method of 
ISarazin fc Tokovinin! I|2002r ). that does not consider the ver- 
tical distribution of the wind speed but just the wind speed 
in two precise regions of the atmosphere (near the ground 
and at 200 hPa), could be as great as 20%-50%, even us- 
ing the appropriated constant. In other words, they proved 
that the proportionality between Vo and V200 is poorly re- 
liable if one wishes precise estimates o f tq (even if we select 
an appropriated constant jj. Masciadri fc Egnerl l|2006l ) con- 
cluded therefore that, for the calculation of the to, the verti- 
cal distribution of the wind speed on the whole troposphere 
is fundamental and necessa ry and the method suggested by 
ISarazin fc Tokovinin! (|2002f ) can provide only some qualita- 
tive (but not accurate) estimates because it presents some 
intrinsic weak points. It appears therefore absolutely impor- 
tant in the field of the site characterization of an astronom- 
ical site, at least in cases such as the calculation of to, to 
characterize the wind speed vertical distribution. How to do 
that? 

The estimate of the wind speed vertical distribution up 
to 20 km at an astronomical site (usually placed on the 
top of mountains) is not trivial. The Generalized Scidar, 
an optical instrument based on a remote sensing princi- 
ple, c an measure the wind speed (jAvila. Vernin fc Sanchez! 
l200ll ) at all the heights in which turbulent layers are present 
i.e. it can reconstruct a sort of vertical profile. Some es- 
timates have been done in the past above different sites 



1 More recently, other authors iGarci'a-Lorenzo et a,l.ll2008f) cal- 
culated Vo = rnax{V groun< i, 0.4 * V200) above the Teide Obser- 
vatory and they found a third different val ue of the constant. 
However, we think they misunderstood the Masciadri & Egncr 
( 
( 



2006) thesis because they stated that the Sarazin fc Tokovinin! 



2002) method needs to be calibrated to be used. 



JAvila et al.ll200d ; lEgner et alj|2007l ; iMasciadri etTaflbOlCh • 
However, this method requires the employment of an in- 
strument that has to be placed at the focus of a telescope 
with a pupil size of at least 1.5 m and it can be performed 
for short periods (typically some tens of nights) related to 
dedicated experiments. Such an instrument is not suitable 
for routinely monitoring the wind speed or climatological 
studies. Alternatively wind speed profiles are routinely cal- 
culated by the General Circulation Models (GCM) (mainly 
from the ECMWF and NOAA/NCEP) and data-set can be 
retrieved in any site in the world with a horizontal resolution 
of 0.25 degrees (operational analyses) and 2.5 degrees (re- 
analysesj 2 ]. Both have been used for astronomical applica- 
tions. More precisely , the ECMWF analyses have been used 
at mid-latitude sites (IMasciadri fc Garfi as 2001; Egner et al.l 
120071 : IMasciadri et al.l 120101) as we l l as at extreme lat - 
itudes JGeissler fc Masciadril l200d ; lHagelin et all 120081 ). 
NOAA/NCEP re-an alyses (vertical pro files) have been used 
at mid-latitude sites (|Carrasco et al.ll2005l : I Avila et al.ll200d : 
iBounhir et"alll2009h and ERA-reanalyses of the ECMWF 
have been used at extreme latitudes (|Sadibekova et al.l 
120061 ). Model outputs showed good correlations with mea- 
surements in all cases in which this has been calculated. The 
unique problem with analyses and re-analyses is that these 
estimates are representative of the wind speed above the 
astronomical site but not in the the surface and boundary 
layer where the local orographic effects have a major effect 
on the wind speed. Above roughly 1-2 km from the ground, 
the wind speed is almost horizontally homogeneous and the 
vertical distribution is basically the same with respect to a 
horizontal extension of some tens of kilometers. Below this 
height the analyses from the GCMs are less representative 
of the wind speed because their horizontal resolution is too 
low to give an accurate description of the interaction of the 
atmospheric flow with the topography near the surface. 

However, in the surface (lowest few tens of meters) and 
boundary layers (typically the first kilomet er above the sur- 
face) it has been proved |Masciadril 120031 ) that mesoscale 
models, with a horizontal resolution of 1 km, can provide 
much better estimates than what analyses from the GCMs 
do above mid-latitude astronomical sites. In that paper the 
author proved that mesoscale models can provide estimates 
better correlated to measurements than analyses from GCM. 
It has also been proved that mesoscale models are able, 
contrary to the analyses from the GCMs, to discriminate 
the wind speed near the ground between two astronomical 
sites (Paranal and Maidanak) characterized by a median 
wind speeds that differ for 4-5 ms~ . This study lets us 
think that mesoscale models could be a useful tool to recon- 
struct the wind speed vertical profile all along the 20 km 
from the ground. They are supposed to be comparable in 



2 To obtain re-analyses the GCM models are re-run as to offer 
an equal horizontal and vertical resolution and the same model 
configuration for many years in the past. In this way, the GCM 
outputs referring to the most recent years and the GCM outputs 
referring to older periods are made uniform. In other words, the 
performances of the re-analyses are degraded with respect to op- 
erational analyses to preserve the temporal uniformity. On the 
other hand the data-set is uniform all along decades and this 
is the reason why re-analyses are used mainly for climatological 
studies. 



Wind speed vertical distribution at Mt. Graham 3 



Table 1. Mcso-NH model configuration. In the second column 
the horizontal resolution AX, in the third column the number 
of grid points and in the fourth column the horizontal surface 
covered by the model domain. 

AX (km) Grid Points Surface (km) 



model 1 
model 2 
model 3 



10 

2.5 
0.5 



80 x 80 

64 x 64 

120 x 120 



800 x 800 

160 x 160 

60 x60 



performances to the ECMWF models above 1-2 km from 
the ground and they are supposed to provide better per- 
formances of the wind speed in the first 1-2 km from the 
ground. 

In this paper we provide a climatological characteriza- 
tion of the wind speed above 1 km from the ground on the 
time scale of ten years using analyses from the ECMWF 
general circulation models and we investigate the possibil- 
ity to use a mesoscale model to systematically reconstruct a 
complete wind speed vertical profile extended on the whole 
20 km above Mt. Graham (Arizona, US) the site of the Large 
Binocular Telescope (LBT). 

In Section 2 we use the operational ECMWF analyses 
over 10 years, 1998-2007, to present a climatological esti- 
mate of the median monthly wind speed. The operational 
data was chosen because of their higher resolution (0.25 de- 
grees) with respect to the re-analyses. We will first prove 
the reliability of the ECMWF analyses in these regions of 
the world comparing analyses with radiosoundings launched 
in from Tucson International Airport ( ~120 km from Mt. 
Graham) . 

In Section 3 we investigate the reliability of the wind 
speed vertical profiles retrieved from simulations with a 
mesoscale model (Meso-NH) in the high as well as in the 
low part of the atmosphere. The mesoscale model is run in 
a grid-nesting configuration covering a total surface of 800 
km x 800 km and three imbricated models with 10, 2.5 and 
0.5 km horizontal resolution (Table[lJ. The innermost model 
covering a surface of 60 km x 60 km. To estimate the model 
reliability we use measurements of the wind speed done by a 
Generalized Scidar run at the focus of the Vatican Advanced 
Technology Telescope (VATT) as well as by an anemome- 
ter located on the roof of the same telescope. This section 
aims to evaluate the possibility to use a mesoscale model 
to characterize the vertical distribution of the wind speed 
extending from the ground up to 20 km above astronomical 
sites. This should certainly represent an extremely valuable 
tool to provide an exhaustive monitoring of the wind speed 
above an astronomical site. 

Finally, in Section 4, we present the conclusion of this 
study. 



2 MONTHLY MEDIAN WIND SPEED 

The ECMWF analyses used in this study are the opera- 
tional analyses, downloaded from the MARS archive|j For a 
more detailed description of the data-set we refer the reader 



to l|Masciadri fc Garfiasl I200H ; iGeissler fc MasciadrH l2006h . 
As explained in the Introduction we only consider the wind 
speed above 1 km from the ground in this section. 

In a previous study (|Egner et al.ll2007l ) done at Mt. Gra- 
ham, it has been proved that independent measurements of 
the wind speed done at the summit of the mountain with a 
Generalized Scidar have a good and small relative discrep- 
ancy of the order of 23% with respect to the ECMWF anal- 
yses extracted from the nearest grid point. Si milar results 
have b een obtained above San Pedro Martir bv lAvila et al.l 
(|200rj) . 

We try here a deeper analysis. The quality of the analy- 
ses in a region of the earth depends on how far/close are the 
meteorological stations with respect to the site one is study- 
ing and on the density of the measurements, i.e the number 
of the meteorological station in that particular region. In 
other words, it depends on the observations that have the 
most important weight in the data assimilation procesQ of 
the GCMs for that region. The nearer to the astronomi- 
cal site are the meteorological stations from which the ra- 
diosoundings are launched, the better is the quality of the 
analyses in proximity to the astronomical site. In order to 
make an estimate of how well the ECMWF analyses perform 
in the region of Mt. Graham we have made a comparison be- 
tween the radiosoundings of the closest available meteoro- 
logical station: Tucson airport (32.23°N, 110.96° W) at ~120 
km southwest from Mt. Graham (32.70°N, 109. 89° W) (Fig[Q- 
top left panel) and ECMWF analyses extracted from the 
closest grid point to the peak of Mt. Graham i.e. (32.75°N, 
110.00°W) at ~ 12 km north-west from the Mt. Graham 
peak (Fig[T]-right panel, point with label (A)). The ra- 
diosoundings are available at 00:00 and 12:00 UTC. We con- 
sidered those at 12:00 UTC corresponding to night time con- 
ditions at local time (12:00 UTC equals 05:00 mst). Fig(2] 
shows the median wind speed vertical profile related to the 
whole year 2005 (left-panel), the difference of the median 
values (central-panel) and the statistical uncertainty a/y/N 
that measures the accuracy of the measurements (right- 
panel). We observe that the difference of the median value 
is always very small (of the order of 1-2 ms _1 ) in most part 
of the 20 km with a relative discrepancy less than or equal 
to 13%. Only at 17 km the difference of the median values 
is somehow larger (~ 4 ms _1 ). However, as can been seen 
in Fig[2]-left panel, at this height the wind speed strength 
is almost 1/3 of the value assumed at the jet-stream level 
and it affects the integrated astroclimatic parameters in a 
less important way. In the calculation of the statistical un- 
certainty a/\ // ~N, N = 328 is the number of nights with ra- 
diosoundings available during 2005. We eliminated 37 nights 
for which the radiosoundings did not cover the whole 20 km 
above sea level. We conclude, therefore, that, above 1 km 
from the ground, the ECMWF analyses provide an accurate 
estimate of the wind speed and the wind speed vertical dis- 
tribution is uniform on a horizontal scale of a few tens of 
kilometers. 

Fig. [3] shows the monthly median vertical wind speed 
profile from the ECMWF-analyses for every month calcu- 



4 The data assimilation is the procedure that, in a forecast model, 
prevents it from drifting away from the true state of the atmo- 
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Figure 1. Top Left: Digital elevation map (DEM) of the Mt. Graham region (model 1 of Table [T] extended on 800 km x 800 km with 
AX=10 km). The distance between the Mt. Graham summit and the Tucson Airport Station (label RS) is ~ 120 km. The black square 
shows the location of model 2. Top Right: DEM (model 2 of Table Q] extended on 160 km x 160 km with AX=2.5 km). The black 
square shows the location of the innermost model (model 3 of Table [l] extended on 60 km x 60 km with AX=0.5 km). The location of 
the Mt. Graham Observatory is marked with an x on the summit of the mountain in the middle of the figure. The grid-points in which 
ECMWF analyses have been extracted (when the GCM resolution is 0.25°) are also marked with a x. The grid-points in which ECMWF 
analyses are extracted (when the GCM resolution is 0.5°) are marked with black dots. Bottom: DEM (model 3 of Table [l] extended on 
60 km x 60 km with AX=0.5 km). 



lated on a time scale of ten years, from 1998 to 2007 at 06 
UTC (23 mst), from the ground up to 25 km. Considering the 
difference in altitude between the ECMWF grid-point (~ 
1320 m) and the Mt. Graham summit (3200 m) we discuss 
these results for h > 4 km above the sea level i.e. roughly 1 
km above the ground of the highest location (Mt. Graham). 
The wind speed profiles calculated above Mt. Graham are 
characteristic of a mid-latitude site and the maximum of 
the wind speed is observed at 11-12 km above sea level. The 
wind speed maximum is located at the height of the jet- 
stream for the most part of the year. During July and Au- 
gust the wind speed at the jet-stream level is much weaker 
than during the other months and, in this period, the high- 
est wind speed value is observed above 20 km, well into the 
stratosphere. The wind speed at the jet-stream level follows 
the classical sinusoidal seasonal trend in different periods of 
the year. Table [2] reports, for each month of the year, the 
average of the wind speed, the standard deviation for the 
yearly average values and the height corresponding to 200 
hPa. 

The wind speed at this height is characterized by values 



comp arable to those observ ed above the major astronomical 
sites (jCarrasco et al.l 120051 ). The month with the strongest 
wind speed at 200 hPa is February (37.21 ms _1 ), while the 
weakest wind speed is observed in July (11.17 ms _1 ). The 
mean wind speed values as well as the standard deviation 
in each month appear very similar to that observed at San 
Pedro Martio This is not surprising considering that the 
two sites are located close to each other. 

The monthly median wind speed profile for each year 
is shown in Fig. [3] The difference between the individual 
years is smallest during the summer, when the strength of 
the wind speed is also the weakest, as is also indicated by the 
smaller difference between the first and the third quartiles 
observed for these months in Fig. [3] The greatest differ- 
ence between individual years occurs during autumn (Octo- 
ber and November) and spring (March and April). The dif- 



5 We precise, for correctness, that studies in ICarrasco et al.l 
(2005) are done with re-analyses while our study is done with 
operational analyses and that we do not study exactly the same 
time period. 
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Figure 2. Left: median wind speed extended on the whole year 2005 calculated for the radiosoundings and the ECMWF analyses (see 
text). Centre: difference of the median wind speed (radiosoundings minus ECMWF analyses). Right: statistical uncertainty calculated 
for N=328, where N is the number of the nights for which we consider radiosoundings. 



Table 2. The mean wind speed at 200 hPa (200 mbar) at Mt 
Graham for the period 1998-2007. Column three and four are the 
standard deviation of the wind speed and the height correspond- 
ing to 200 hPa. 



Month 


Avg (ms 1 ) 


std (ms 1 ) 


h (km) 


January 


32.79 


4.22 


11.92 


February 


37.21 


4.55 


11.88 


Mars 


32.65 


6.70 


11.89 


April 


36.09 


6.81 


11.98 


May 


26.77 


6.54 


12.11 


June 


21.64 


5.36 


12.26 


July 


11.17 


2.08 


12.37 


August 


11.80 


2.82 


12.37 


September 


22.77 


3.42 


12.31 


October 


29.79 


5.85 


12.16 


November 


30.88 


6.25 


12.04 


December 


32.21 


3.20 


11.93 


Average 


27.06 


4.82 


12.10 



ference between the calmest and the strongest yearly wind 
speed at the jet-stream level is over 20 ms -1 during these 
months. This is coherent with the dispersion (dashed lines) 
indicated in Fig. [3] A peculiar pattern of alternating years 
with stronger and weaker winds is also found in October. 
The other months do not show this pattern, but there are 
clearly rather large differences in-between different years. 
This result indicates that studies on seasonal trends, for the 
wind speed and also the optical turbulence (that depends 
on the wind speed), should be preferably done on time scale 
of the order of some years to filter out effects due to this in- 
trinsic variability that characterizes the wind speed in each 
month. 

It is worth noting that the horizontal resolution of 
ECMWF-analyses, has changed during the studied period 



(1998-2007]j. During this period the horizontal resolution 
changed from 0.5 ° to 0.25°. Figure [l] shows the locations 
of the four grid points that are closest to Mt. Graham in 
the case of a resolution of 0.25° (crosses) and 0.5° (black 
points). With a horizontal resolution of 0.25°, the closest 
grid point to the peak of the mountain is A=(32.75°N, 
110.00°W), located west-northwest of the peak. With a 
horizontal resolution of 0.5° the closest grid point is B= 
(32.50°N, 110.00°W), southwest of the peak. A compari- 
son of the vertical wind speed profile from both these grid 
points (see Appendix A) show that they are almost identi- 
cal in the free atmosphere, therefore we consider that there 
are no problems of data inhomogeneity when we treat the 
average on time scale of the order of ten years. 



3 THE WIND SPEED FROM THE MESO-NH 
MODEL 

The Meso-NH is a non-hydrostatic mesoscale model devel- 
oped jointly by Me teo-France and Laboratoire d'Aerologie 
(Lafo re et al.lll998l ). It is a grid point model based on the 
anelastic approximation that can simulate the temporal evo- 
lution in three dimensions of the classic meteorological pa- 
rameters such as wind speed and direction, potential tem- 
perature and pressure. 

For this study we have run the Meso-NH model in grid- 
nesting mode, using three two-way nested models, centered 
at the peak of Mt Graham. Details of the horizontal size 
and resolution of the three models are reported in Table [T] 
The outermost model covers an area of 800 x 800 km. The 
areas covered by the three models are shown in Fig. [T] The 
vertical grid is composed of 49 levels covering up to 20 km 
above sea level. The first vertical grid point is located at 20 
m. Above we have a logarithmic stretching of 20% for the 



http : //www . ecmwf . lnt/products/data/technlcal/model_id/lndex . html 
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Figure 3. The median wind speed calculated with the ECMWF analyses for each month during 1998-2007 at Mt Graham, solid line. 
The dashed lines represent the first and third quartile. The dotted line indicates 200 hPa. Data-set are extracted from the grid point 
indicated with the label (A) in Fig[l}right panel. 



vertical grid size up to 3500 m_|. Above 3500 m the distance 
between the vertical levels is fixed and equal to 600 m. 

The model is initialized with the analyses of the 
ECMWF, which also provides the boundary conditions for 
the outermost model. The runs start at 00:00 UTC (17:00 
mst), the synoptic hour closest to local evening, and last 
for 12 hours, to early morning in Arizona (05:00 mst). The 
first two hours are rejected to avoid the calculations be- 
ing affected by spurious values due to the adaptation of the 
atmospheric flow to the ground. Data from the innermost 
model with the highest horizontal resolution (model 3) are 
treated here. 

To ensure that the Meso-NH model can reconstruct re- 
liable wind speed profile for h > 1 km we compare the wind 
speed calculated by Meso-Nh with the wind speed calcu- 
lated by the ECMWF analyses in the four grid points that 



7 The logarithmic stretching creates a vertical grid where the dis- 
tance between each level is 20% more than between the previous 
two. 



surround Mt. Graham using the 0.25° resolution (the four 
grid points are indicated with a £ in Fig. [T]). Figure [5] shows 
this comparison calculated for one night, May 21, 2005. The 
Meso-NH profiles are obtained by averaging the outputs cal- 
culated at each hour for the entire 10 hour simulation cover- 
ing the period (17:00 to 05:00 mst). The ECMWF data are 
only available at the synoptic hours, therefore we considered 
the average of the ECMWF data at 00, 06 and 12 UTC (at 
17, 23 and 05 mst). 

We observe that the wind speed profile from Meso-NH 
is very well correlated with the ECMWF analyses for all the 
four grid points above roughly 1 km from the ground. This 
also indicates that the wind speed is horizontally homoge- 
neous over a surface of 0.25° x0.25°. 

What is the Meso-NH ability in reconstructing the tem- 
poral evolution of the wind speed? To answer to this ques- 
tion we choose two nights, one characterized by a weak wind 
speed (May 29, 2007) and another characterized by a strong 
wind speed (February 26, 2008), and we compared the wind 
speed profiles calculated with the ECMWF analyses and the 
Meso-Nh model at the start and at the end of the Simula- 
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Figure 4. The temporal evolution of the annual median wind speed on the 1998-2007 period calculated with ECMWF analyses at Mt 
Graham. The color scale goes from ms" 1 (blue) to 40 ms" 1 (dark red). 



tion. We consider the wind speed calculated by the Meso- 
NH at the peak of Mt. Graham and the ECMWF wind speed 
extracted in the nearest grid point to the Mt. Graham sum- 
mit. Again we are only interested, in this phase, in the wind 
speed above 1 km from the ground and, following the same 
logic we used in Section [2] we discuss the results obtained 
above 4 km from the sea level. Figure [6] shows the result 
of the comparison. ECMWF profiles are interpolated to the 
Meso-NH vertical grid points. Fig. [6}top and Fig. [6]-bottom 
show, respectively, the start and the end of the simulation 
We observe that the wind speed profile calculated by Meso- 
Nh has evolved during both of the nights accordingly with 
the wind speed as calculated by the ECMWF. The shape 
of the profile is reasonably well correlated as well as the 
strength of the wind speed. We observe that the wind speed 
increases somewhat during the night on May 29, 2005 while 
decreases its strength on February 26, 2008. We conclude 



therefore that the mesoscale model is able to reconstruct 
the wind speed vertical distribution in the high part of the 
atmosphere and reproduces the spatio-temporal wind speed 
variability in a satisfactory way. In the case of May 29, 2007 
(Fig. [BJbottom) the wind speed in the low part of the at- 
mosphere reconstructed by the Meso-Nh model presents a 
few peaks that might be originated by the better horizontal 
resolution of the mesoscale model. 

What about the ability of Meso-Nh in reconstructing 
the wind speed in the first kilometer from the ground? For 
this study we consider measurements provided by two differ- 
ent instruments: (1) measuremen ts of the wind speed p rofile 
made with a Generalized Scidar (Masciadri ct al. 2010) and 
related to 39 nights in different seasons (Table[3]). (2) 'in situ' 
measurements of the wind speed done with an anemometer 
located on the roof of the VATT (~ 20 m from the ground). 

For these nights we can retrieve from the GS the wind 
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Figure 5. The vertical wind speed profile for the four grid points from the Meso-NH (solid line) and the ECMWF analyses (dashed 
line) that are located closest to Mt. Graham (grid points indicated with a cross in Fig[TJ, when using the 0.25°resolution. Calculations 
are done for the night May 21, 2005. The dotted line indicate the typical height below which the comparison makes no sense. 



speed along the whole 20 km. As already mentioned, a com- 
parison of the Generalized Scidar wind speed profiles and 
the ECMWF ana lyses at Mt, Graham has been done in 
lEgner et al.l (|2007l ) and it has been found a good correla- 
tion for h > 1 km. Qualitative comparisons of the Gener- 
alized Scidar wind profiles with NCEP/NCAR re-analysis 
and wind speed from in situ balloons for 15 ni ghts in May 
2000 at San Pedro Martir has been described inlAvila et al.l 
IJ2006h . Also iGarcia- Lorenzo fc Fuensalidal l|200fj ) compared 
the Generalized Scidar wind speed measurements done at 
the Teide Observatory with radiosoundings data launched 
by a meteorological station placed at ~ 13 km away from 
the summit for four summer nights in 2003. However, to our 
knowledge, there is no detailed study of the reliability of the 
Generalized Scidar at describing the wind speed in the sur- 
face or boundary layer. For this reason, the measurements 
from an anemometer have been used to have an independent 
measurements in the very low atmosphere. 

We have calculated the average wind speed vertical pro- 
file in the first kilometer for all 39 nights for which there are 
measurements of the wind speed from the Generalized Sci- 
dar (Fig. [7] - blue line) and compare the resulting profile 
with the average wind speed simulated by the Meso-NH for 
the same nights (Fig. [7] - green line) . We remind that the 
GS is placed at the focus of the VATTjfl therefore the first 



8 The Vatican Advanced Technology Telescope is also located at 
the top of Mt. Graham at around 250 m from the LBT. 



Table 3. Observing runs at Mt. Graham and the number of nights 
for which there are wind speed profiles from the GS. 



Obs. runs (utc) 



Nights 



27 April 2005 


1 


20-26 May 2005 


6 


7-16 December 2005 


5 


28 May - 4 June 2007 


8 


17-29 October 2007 


11 


24 February - 4 March 2008 


8 



measurement (on the vertical grid) for the GS is at ~ 20 
m from the ground. Looking at Fig[7]we observe that, al- 
most everywhere the two profiles are very well correlated, 
with the Generalized Scidar estimating wind speed that dif- 
fers from the Meso-Nh model for not more than 1 ms _1 all 
along the first kilometer from the ground. At 20 m the GS 
wind is slightly stronger than the wind calculated by the 
model. This is probably due to the fact that the GS mea- 
surements have a vertical error bar of the order of 25-30 m 
(see Masciadri et al. 2010) and the wind shear is particu- 
larly strong at this height. Therefore it might be that the 
higher wind speed detected by the GS at h =20 m can be 
originated by thin layers flowing at slightly higher distance 
from the dome and not resolved by the instrument. Figure 
[8] shows the average wind speed near the ground for every 
night reported in Table|3]in chronological order as measured 
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Figure 6. Vertical wind speed profile of the Meso-NH (solid line) and the ECMWF analyses (dashed line) at the beginning (00:00 UTC) 
and the end (12:00 UTC) of a simulation for two nights, May 29, 2007 (left hand side) and Feb. 26, 2008 (right hand side). 
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Figure 7. Left: Mean vertical profile of the wind speed in the boundary layer (1 km) obtained with the GS (blue line) and the Meso-NH 
(green line). The error-bars show the standard deviation. Right: zoom in the first 200 m. 



by the Generalized Scidar and reconstructed by Meso-NH at 
around 20 m i.e. at the height of the dome of the telescope. 
In the same figure is also reported the wind speed measured 
in the same nights by the anemometer mounted on the top 
of the VATT. From a qualitative point of view, in Fig. [S] 
we observe that the temporal evolution of the wind speed 
reconstructed by the model during the 39 nights is very well 
correlated to the anemometer trend. The wind speed recon- 
structed by the model decreases and increases following the 
wind speed evolution measured by the anemometer in the 
same nights. Only a small off-set is present between the two 



estimates. From a quantitative point of view the mean wind 
speed from the GS at h = 20 m is 7.37 ms _1 , the mean wind 
speed reconstructed by the Meso-Nh model is 5.21 ms _1 , the 
mean wind speed measured by the anemometer is 2.44 ms~ . 
The wind speed reconstructed by the model is well included 
in the range of the wind speed measurements and this cer- 
tainly proves the reliability of the calculated wind speed. 
The weaker wind speed from the anemometer is probably 
due to the fact that the anemometer is an 'in situ' measure- 
ment that is done well below the top of the trees. It has been 
observed in the past that the friction of the atmospheric flow 
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Figure 8. The average wind speed near the surface for Meso-NH 
(thick solid line and triangles) at h=20 m, an anemometer placed 
on the roof of the VATT at h ~ 20m (dashed line and stars) and 
the Generalized Scidar (thin solid line and circles). 



with the trees causes a sharp and drastic decreasing of the 
wind speed below this height (something that is confirmed 
in the profile reconstructed by the model below 20 m). 

We conclude therefore that the model provides reliable 
estimates all along the whole 20 km. In the future it would 
be good to be able to equip the observatory with anemome- 
ters located at different heights below and above the top 
of the trees preferably in open space environment so that 
the wind speed is not affected by the presence of buildings. 
This should permit us to better monitor the particular sharp 
change of wind speed in this region. 



4 CONCLUSIONS 

In this paper we aimed to give a complete characterization 
of the wind speed vertical distribution at Mt. Graham (Ari- 
zona, US) for astronomical applications. The simplest way 
of retrieving a complete characterization of the wind speed 
profile (with exception of the surface and boundary layer) is 
using the data from General Circulation Models (GCMs). In 
this study we use the operational analyses from the ECMWF 
extracted from the grid point closest to the peak of Mt. Gra- 
ham, with a 0.25 degrees resolution, to study the vertical 
wind speed distribution over 10 years (1998-2007). 

We have verified that the wind speed profile retrieved 
from the operational analyses of the ECMWF model is con- 
sistent with what is obtained by the radiosoundings from the 
nearby Tucson International airport. We also proved that 
the wind speed in this region is homogeneous with respect 
to horizontal spatial scales of the order of some tens of kilo- 
meters above ~4 km from the sea level. 

Having proved that the ECMWF analyses are reliable, 
we presented the monthly median wind speed extended on 
a 10 years time scale. The wind speed profiles are rather 
typical of a mid-latitude site with a pronounced wind speed 
maximum at the jet stream level (10-12 km above sea level) 
during most of the year. On the contrary, during the summer 
the maximum wind speed is located well into the strato- 
sphere. The strongest variability (from different years) in 



the monthly median wind speed is found during spring and 
autumn. 

For the same period we also provided the monthly mean 
wind speed values at 200 hPa, corresponding to the maxi- 
mum wind speed values at the jet-stream height. The month 
with the strongest wind speed at 200 hPa is February (37.21 
ms _1 ), while the weakest wind speed is observed in July 
(11.17 ms _1 ). Results indicate that the wind speed at the 
jet-stream level is very similar to what has been observed 
above the Observatory in San Pedro Martir (Baja Califor- 
nia) and it is consistent with the values observed above the 
best astronomical sites in the world. 

Besides, we proved the reliability of a mesoscale model 
(Meso-Nh) in reconstructing the wind speed on the whole 
20 km just above the summit of Mt. Graham included the 
boundary layer and the surface layer. We proved that the 
wind speed reconstructed by the model is very well corre- 
lated to the ECMWF analyses and it is also able to provide 
realistic profiles in the first kilometer from the ground. 

The wind speed estimates from the mesoscale model 
have been compared to measurements from a Generalized 
Scidar and an anemometer located at 20 m from the ground 
on a sample of 39 nights. Above 50 m the wind speed profiles 
reconstructed by the model match in a very satisfactory way 
(AV < 1 ms _1 ) with respect to the measured wind speed 
profiles. Closer to the surface, just in proximity of the top 
of the trees, the wind speed estimated by the model is in- 
cluded in the range of values given by the anemometer and 
the Generalized Scidar and for this reason can be considered 
satisfactory. However, the dispersion between the anemome- 
ter and the GS measurements seems a little too large. This 
difference is highly probably due to the fact the anemome- 
ter measures a wind speed that is not completely in free air 
(the anemometer is placed beside the dome of the VATT) 
and below the top of the trees. On the other side, the GS is 
probably affected by the wind speed just above the top of 
the trees because of its finite vertical resolution. The quali- 
tative trend of the wind speed observed all along the sample 
of 39 nights is however very well reconstructed by the model 
and in agreement with measurements. The model appears 
to reconstruct very well the wind speed behavior above and 
below 20 m. 

This paper therefore validates the Meso-Nh model as 
a tool to predict the wind speed vertical profile V(h) from 
the ground up to 20 km above Mt. Graham for each night 
and, at present time, it appears as the unique method to 
systematically estimate the whole wind speed vertical profile 
above an astronomical Observatory. We remind that it has 
been proved (Masciadri et al. 1999a, Masciadri & Jabouille, 
2001, Masciadri et al. 2004, Lascaux et al. 2010) that Meso- 
Nh can provide reliable C% profiles above an astronomical 
site and it appears therefore as an extremely useful tool 
for tq estimates. Besides there are evidences that it would 
be very useful to supply the Mt. Graham Observatory with 
anemometers located at different heights below and above 
the top of the trees (~ 20 m) because this should permit to 
provide a better constraints of the model for dedicated and 
more detailed applications. 
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APPENDIX A: THE WIND SPEED AT 
DIFFERENT HORIZONTAL RESOLUTION 

The horizontal resolution of the operational analyses of the 
ECMWF has changed during the ten year-period studied in 
this paper. This implies that the position of the closest grid 
points to Mt. Graham changes. The data in this study are 
downloaded from the 0.25° resolution where the closest grid 
point to Mt. Graham is located east northeast of the moun- 
tain peak (Fig. [[]- point (A)) at 32.75°N, 110.00° W. Using 
the 0.5° resolution the closest grid point is located 32.50°N, 
110.00° W, southwest of the mountain (Fig.[T]- point (B)). 

To closer examine which impact the different resolutions 
have on the vertical wind speed profile we have downloaded 
the data from both of these grid points for the entire year 
2002. The monthly median wind speed profile is presented in 
Fig. IA1I where the data from the grid point using the higher 
resolution (32.75°N, 110.00°W) is plotted with a solid line 
and the data from the 0.5° resolution (32.50°N, 110.00° W) is 
plotted using a dashed line. The difference between the two 
data-sets is very small. During all months the lines over- 
lap each other almost entirely. Some smaller offsets exist, 
but they are generally minor. The largest difference found, 
near the jet stream-level in October, is 2.8 ms~ . We con- 
clude therefore that the change in horizontal resolution in 
ECMWF-analyses during the 10 years did not introduce any 
biases in our calculation. 
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Figure Al. The median monthly wind speed for 2002 at the two different grid points closest to the peak of the mountain using 0.25° 
resolution (32.75°N, 110.00°W) (solid lines) and 0.5° resolution (32.50°N, 110.00°W) (dashed lines). 



